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Abstract. We study modified theories of gravity of the f(R) type in Palatini formalism. We first consider the stability of 
atoms when the Palatini gravitational interaction is taken into account in the derivation of the non-relativistic Schrodinger 
equation. We show that theories with infrared curvature corrections are ruled out by the mere existence of atoms. In particular, 
we carry out fully perturbative calculations that, for the first time, convincingly rule out the 1 /R model of Carroll et al. in 
its Palatini version. We then study the Planck scale corrected quadratic model f(R) = R + R 2 /Rp and show that it can avoid 
the big bang singularity for matter sources which satisfy all the energy conditions. We comment on the mechanisms that cure 
this singularity and point out that they are closely related to non-perturbative terms also present in the atomic Hamiltonian of 
infrared corrected models. 

Keywords: Modified Gravity, Bouncing Cosmology, Backreaction 
PACS: 04.50.Kd, 98.80.-k 

INTRODUCTION 

The currently accepted cosmological model assumes that the matter and energy contents of the Universe are vastly 
dominated by two dark sources whose presence can only be inferred from their gravitational interaction [1, 2, 3]. 
And this gravitational interaction is dictated by Einstein's equations. Obviously, if the dynamics were different, the 
inferred amounts of dark matter and dark energy would be different from those predicted by Einstein's theory. This 
simple observation has led many authors to consider modifications of the gravitational field equations which could 
account for the cosmic speed up without the need for sources of dark energy . One of such approaches assumes that 
the gravity Lagrangian could be some to-be-determined function f(R) which departs from General Relativity (GR), 
f(R) = R, in certain regimes (see [4, 5] and references therein). An interesting aspect of considering modified theories 
of gravity with non-linear curvature corrections in the Lagrangian is that the derivation of the field equations requires 
to specify whether the connection is dependent or independent of the metric. The latter case is known as metric-affine 
formalism, and leads to field equations inequivalent to those found in the more familiar metric formalism, in which 
the connection is assumed to be compatible with the metric, i.e., it is represented by the Christoffel symbols. When 
in a metric-affine theory the independent connection is not coupled to the matter fields, then the theory is said to be 
formulated in the Palatini formalism. 

The purpose of this work is twofold. First we want to show that Palatini f(R) models with correcting terms that 
dominate at low curvatures are ruled out by the very existence of atoms. We will see that the characteristic cosmic 
density scale of such models, which sets the scale at which cosmic speedup arises, can be effectively reached in 
microscopic systems such as atoms. At such density scales, the modified gravitational dynamics induces strong 
backreaction effects in the metric which affect the atomic Hamiltonian and end up disintegrating the atom. This will 
lead us to the second point. The strong backreaction effects that spoil infrared corrected models turn into a very 
effective mechanism to remove singularities in ultraviolet modified models. In fact, we will show that in simple 
models such as f(R) — R + R 2 /Rp, where Rp ~ lp 2 is of the order of the Planck curvature, the Big Bang singularity 
of General Relativity can be replaced by a bounce even if the universe is dominated by matter sources which sat- 
isfy all the energy conditions. If this theory is formulated in the usual metric formalism, no bouncing solutions exist [6]. 

The content is organized as follows. We first derive the field equations of a generic f(R) Lagrangian and briefly 
comment on its weak field limit. Second, we discuss the effective Hamiltonian that governs the motion of an electron 
in a Hydrogen atom and which is subject to the Palatini gravitational interaction. We show that infrared curvature 
corrections in the Lagrangian generate instabilities which end up disintegrating the atom, which clearly rules out such 



models. We then show that in ultraviolet corrected models similar mechanisms also appear in critical scenarios such 
as the very early Universe. The effect in this case is the avoidance of the Big Bang singularity by the emergence of a 
cosmic bounce due to strong backreaction effects. In these models atoms are as stable as in GR. We end with a brief 
summary and conclusions. 



FIELD EQUATIONS 

Let us begin by defining the action of Palatini theories 

S\g,T, Wm] = ^2 / d 4 xV=gf(R) +S m [8nv,Ym] (D 

Here f(R) is a function of R = g^R^T), with fy v (T) given by fy v (T) = -d„rjj; v + d x T* v + 1*,^ - 1*,^ 
where is the connection. The matter action S m depends on the matter fields y m , the metric 8nv, which defines the 
line element ds 2 = g fiv dx fl dx v , and its first derivatives (Christoffel symbols). The matter action does not depend on 
the connection T^ v , which is seen as an independent field appearing only in the gravitational action (this condition 
is not essential and can be relaxed at the cost of introducing a non-vanishing torsion). Varying (1) with respect to the 
metric g^ v we obtain 

f R (R)R^(r)- 1 -f(R)g^ = K% v (2) 

where fg (R) = d f/dR. From this equation we see that the scalar R can be solved as an algebraic function of the trace 
T. This follows from the trace of (2) 

f R (R)R-2f(R) = K 2 T, (3) 

The solution to this algebraic equation will be denoted by R = &(T). The variation of (1) with respect to leads to 
[V~8fn8^ v ] — 0, where = /r(^[T]) is also a function of the matter terms. This equation leads to 

r Aiv = — (<V P v + d v t pfl - dpt^y) (4) 

where t^ v = fygpv, and = f R \^\o^ is dimensionless and normalized to unity outside of the sources (T — 0). It is now 
useful to rewrite (2) adding and subtracting Jf&(T)g^ v = ^ft a ^R a p (r)f jUV to get 

fRG^iv(t) = K 2 Tp V - - — ^ — t^v (5) 

where G^ v (t) is the Einstein tensor associated to t^ v . The equations of motion (5) for the auxiliary metric t^ v are 
(formally) considerably simpler than those for g^: 
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Solving for t^ v using the system (5) and then going back to g^ v via the conformal transformation g^ v = ^)(T)^ l t^ v 
is a useful simplification that makes the task of finding solutions much easier. Note that this fortunate circumstance 
is due to the fact that the conformal transformation completely cancels out the disturbing derivatives on the right 
hand side of (6). Moreover, the conformal relation between the two metrics puts forward the fact that the metric g^ v 
receives two kinds of contributions: the usual non-local contributions that result from integration over the sources 
(like in GR), which produce the term t^ v , and local contributions due to <j)(T), which depend on the local details of 
T at each space-time point. This local contribution arises due to the independent character of the connection and, to 
our knowledge, is not present in any other metric theory of gravity, where the connection is generally assumed to be 
metric compatible (Levi-Civita connection). 

The similarity between the field equations of GR and (5) suggests that for weak sources and reasonable choices of 
Lagrangian f(R) (those that lead to negligible cosmological term ^2/~^ )» the right hand side of (5) is small and, like 



in GR, fuv can be expressed as t liv {x) w r]^ v + h^ v (x), where |/i^ v (^)| <C 1 is given as an integral over the sources 
(for details of exact calculations see [13]). For very weak sources such as atoms or elementary particles, we find that 
the self-contribution to h^y — > 0. On the other hand, if our microscopic system is placed in an external gravitational 
field whose contribution to t^ v is not negligible, we can always take a coordinate system in which the metric at the 
boundaries of a box containing the system (large as compared to the microscopic system but small as compared to the 
range of variation of the external metric t^ v ) becomes <~ 77^ v . In both situations, the metric g^y becomes simply 

gfivW-^^l^V (7) 

where <j>(T) — ► 1 in the (empty) boundaries of our auxiliary box but depart from unity in the region supporting the 
sources (which are assumed to be fields spread over the space). For a detailed discussion of the Newtonian and 
post-Newtonian limits of f(R) theories see [7, 8]. 



ATOMIC INSTABILITIES 

It is well known [9, 10, 11] that the energy levels of a Hydrogen atom falling freely in an external gravitational field (in 
GR) will be shifted in a very characteristic way due to the interaction of the electron with the curvature of the space- 
time. Though external fields in Palatini theories of gravity must also lead to this phenomenon, we will focus here on a 
different aspect. We will study the effect that the local energy-momentum densities have on the non-relativistic limit 
of the Dirac equation due to the factor (T) ~ 1 appearing in (7). The non-relativistic Schrodinger equation that governs 
the motion of an electron in an external electromagnetic field taking into account the Palatini gravitational interaction 
is the following (see [12, 13] for a detailed derivation and discussion) 

gr\ = 1^1 — [(p-eA) 2 -ea-B}+eA \ri 
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where $ represents the total energy minus the rest mass energy mo, Q. = (3/4)ln0(r), m = m<j>~?, m <~ mo is the 
mass appearing in the action of the Dirac spinor, and 77 is the large component of the Dirac spinor (the positive energy 
Foldy-Wouthuysen bispinor). 

The first line of this equation is very similar to the well-known non-relativistic Schrodinger-Pauli equation (see (10) 
below). The only difference being the term 1 / (m + mo). The second and third lines, however, represent completely 
new terms generated by the Palatini gravitational interaction. When the gravity Lagrangian is that of GR, 0(r) = 1, 
we recover the Schrodinger-Pauli equation if mo is identified with m. 

To extract physical results from this equation, we proceed as follows. We first solve (8) in the case of GR, 
f(R) — R,(j)(T) = 1, which is well known. Then we switch to a different gravity Lagrangian (assuming that we 
have the ability to do that) and study how the system reacts to that change. The reason for this is that in a general 
f(R) the metric is sensitive to the local T^ v via 0(r) _1 , and changes in the metric due to the matter distribution 
could react back on the matter equations. If the new interaction terms in (8) lead to small perturbations, then the 
initial wavefunctions will be, roughly speaking, stable with perhaps small corrections which could be computed using 
standard approximation methods. If, on the contrary, the energy associated to the gravitationally-induced terms is 
large, that would mean that the original configuration is not minimizing the modified Hamiltonian and, therefore, 
large modifications would be necessary to reach a new equilibrium configuration. Depending on the magnitude of the 
reaction on the system, we could estimate whether the theory is ruled out or not. 

4 

Let us now focus on the infrared corrected model f(R) = R — initially proposed in [14] within the metric 
formalism and in [15] in the Palatini version, which is characterized by the low curvature scale ji 2 . In this particular 
model, we find that 3S{T) = -(k 2 T + ^(k 2 T) 2 + 12jU 4 )/ 2 recovers the GR limit for |fc 2 r| > jx 2 and tends to a 



constant 8% ~ jx 2 for | K 2 T\ <C /I 2 . It also follows that the function <j> (T) is given by 
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where x = —T/T c , T c 
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; /I 2 / K" 2 = p^, and <~ 10~ 26 g/cm 3 represents the characteristic cosmic density scale of the 
theory, which triggers the cosmic speedup. It is easy to see that at high densities, as compared to p^, <j)(T) — ► 3/4, 
whereas for p < p M we find <j>(T) — > 1. 

Expressing length units in terms of the Bohr radius («o <~ 0.53 • 10~ 10 m), we find that X = = 10 24 P e (x), where 
we have intentionally omitted the nuclear contribution (only relevant at the origin) for simplicity and have used 
p e {x) = mP e (x) with P e (x) = TJ* (x)r](x). This expression for x indicates that the electron reaches the characteristic 
cosmic density, x <~ 1, in regions where the probability density is near P e (x) ~ 10~ 24 . In ordinary applications, one 
would say that the chance to find an electron in such regions is negligible, that that region is empty. In our case, 
however, that scale defines the transition between the high density (t >• 1) and the low density (t <C 1) regions. 
In regions of high density, we find that rapidly tends to a constant, 0o» = 3/4, which leads to m — 2m/\/3 and 
Vi2 = 0. If we then identify m — ► \/3mo/2, equation (8) reduces to the usual Schrodinger-Pauli equation 
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This fact clarifies the introduction of mo above and guarantees that the solutions of the equation in this region coincide 
with the standard ones (those of GR). Therefore, the perturbative approach that we have followed is well justified. Let 
us now see what happens in regions of low density. In those regions, <j)(T) tends to unity, Vi2 = 0, and m — > m as x — > 0. 
The mass factor dividing the kinetic term is now a bit smaller (mo > m) than in the high density region. But the mass 
difference m — mo is not zero. This is a remarkable point, because m — mo ~ — 0.1 3mo is negative and of order <~ mo, 
which represents a large contribution to the Hamiltonian. Modifications of the standard solutions are thus expected in 
this region. To better understand the effect of this term, it is useful to consider the ground state, J](i o o) — e / (g> 

where \ \,s) represents a normalized constant bispinor. In this case, the transition from the high density region to the 
low density region occurs at r rts 26oq. In Fig.l we have plotted the most representative potentials in dimensionless 
form 
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where V e is the electrostatic potential generated by the proton, lengths are measured in units of the Bohr radius, 

13.6 eV. Note that V& = V 2 i2 — |Vi2| 2 only contains the most important 



x = r/ao, and energies in units of 

contributions associated to Q.. 

In this case, Vn represents a small transient perturbation. The mass difference V m , however, introduces a deep 
potential well in the outermost parts of the atom that must have important consequences for its stability. Note that 
this effect is not an artifact of the non-relativistic approximation, since it also occurs in the full relativistic theory 
(see [12] for details). In the initial configuration of the atom, corresponding to GR, the wavefunction of the ground 
state is concentrated near the origin, where the attractive electric potential is more powerful (V e — ► — °°). As we switch 
on the l/R theory, a deep potential well of magnitude ~ — 0.1 3mo appears in the outer regions of the atom, where 
p e (x) < Pji, which makes the ground state unstable and triggers a flux of probability density (via quantum tunneling) 
to those regions. The half life of Hydrogen subject to this potential can be estimated using time dependent perturbation 
theory (see [12]) yielding 

h 



- « 6 • 10 J 



(14) 



We thus see that the initial, stable configuration is destroyed in a lapse of time much shorter than the age of the 
Universe, which is in clear conflict with experiments. 
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FIGURE 1. Contribution of the different potentials in the ground state. The solid line, which represents the sum of all the 
potentials, tends to the constant value — 0.134m ( ) (or —5048 in the units of the plot). 
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FIGURE 2. The different contributions in this plot are V e ~ -1, V m ~ -5 • 10 3 , and V a ~ ±10 19 . The y-axis is measured in units 
of 13.6eV; the x-axis in units of a - 



Further evidence supporting the instability of the atom is found in the existence of zeros in the atomic wavefunctions 
in between regions of high density because, obviously, before (and after) reaching p e {x) = the characteristic scale 
p e (x) <~ Pfi is crossed. The first excited state, J](2oo) = / 1 (1 ~ "£~) e ^^ 2 " ® \\i s )> nas a zero at r = 2ao- The 

radial derivatives of 0(r) at that point are very large and lead to very important perturbations which overwhelmingly 
dominate over any other contribution (see Fig.2). 

The magnitude of Va — V 2 i2 — |Vi2| 2 at r = 2oq oscillates between 10 20 and — 10 21 eV in an interval of only 
2 • 10~ 10 ao. Needless to say that this configuration cannot be stable and that strong changes must take place in the 
wave function to reduce the energy of the system. Such changes should tend to reduce the magnitude of the density 
gradients (Vi2) to minimize the value of Va, which will likely lead to a rapid transition to the ground state, where 
Vh is small. One can easily verify that strong gradients Vi2 also appear at the zeros of all the tj B) o,o wavefunctions, 
which generate large contributions Vh in those regions. Furthermore, if one considers stationary states with I ^ 0, Va 
has important contributions not only at the zeros of the radial functions, but also at the zeros of the angular terms. 
Thus, the pathological behavior described for the spherically symmetric modes gets worse for the 1^0 states. One 
thus expects the decay of these states into states with less structure (weaker gradients) such as the ground state, which 
will later decay into the continuum. All this indicates that the existence of bound states, with localized regions of 
high probability density (where "high" means above the scale p^), are impossible in this theory because of the large 
gradient contributions Va and the deep potential well V m . 

Though we have only analyzed in detail the infrared-corrected model f(R) =R — jJ. 4 /R, the instabilities associated 
to the potential well m — mo and the zeros of the wavefunction must be present in all gravity models sensitive to 
low curvature/energy-density scales. Therefore, such models are ruled out by the mere existence of atoms. Note that 
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FIGURE 3. Bouncing cosmologies with K = compared with the GR solution (red lines). The continuous black curve represents 
the case w = 0, corresponding to f(R) = R — R?/\Rp\. The dashed line represents the case w = 1 of the theory f(R)=R + R 2 /\R P \. 
The initial density is the same for the four curves. 

models with high curvature corrections, such as f(R) = R+R 2 /Rp, do not exhibit any instabilities of those affecting 
infrared corrected models. In this example, we find that = 1 — 2k 2 T/Rp, which implies that all gradients appearing 
in Vq are suppressed by the Planck density pp = Rp/ K 2 and also that m = mo(l + 0\p/pp}) ~ my. The effects of the 
potentials induced by the (j) (T)~ 1 factor in front of the metric are, therefore, much smaller than the corrections due to 
the Newtonian potential. Thus, in such models atoms are as stable as in GR. 

RESOLUTION OF COSMIC SINGULARITIES 

The unexpected effects of the Palatini gravitational interaction for infrared corrected models in atomic systems raises 
a number of natural questions. What is the physical mechanism responsible for the atomic disintegration? Where does 
that energy come from? Understanding these points is very important because in the atomic system studied above 
the gravitational field has been able to surpass the electrostatic forces that keep the electron bound to the nucleus. 
This means that we are facing a new gravitational mechanism able to generate strong repulsive forces even when tiny 
amounts of matter and energy are involved. It should be noted that such forces arise in regions where the gradients 
d^fp and V^V v /r are strong. In fact, it is precisely those terms which represent the maximum contribution to the 
Hamiltonian through the induced potential Vq and which introduce strong non-perturbative effects in the atomic 
dynamics. Thus, the rapid variation of the terms d^fp and V^Vy/^ in regions where the characteristic scale of the 
Lagrangian is reached generates important contributions on the right hand side of the metric field equations (6) and 
results in a strong backreaction. The question now is to see if such repulsive forces and backreaction effects can be 
achieved in other scenarios of interest. We will show next that ultraviolet corrected models may also give rise to very 
strong repulsive forces in critical situations such as in the very early Universe. This will confirm that Palatini f(R) 
theories possess the necessary ingredients to get rid of the Big Bang singularity in quite general situations. 

As an illustrative example, let us consider the following modified Lagrangian f(R) =R + R 2 /R P , where R P 
represents a curvature scale which could be identified with the inverse Planck length squared. For this Lagrangian we 
find that 0?, = — K 2 T and = 1 — 2k 2 T /Rp. A glance at the field equations (6) indicates that the right hand side is 
dominated at low curvatures by the usual Einstein term K 2 T^ V , with corrections of the order ~ K 2 T /Rp which only 
become important at very high curvatures. At low curvatures/densities, the dynamics is thus very well described by 
GR. It is important to note that the terms d^ fpdvfp are of order ~ (k d^T /Rp) 2 and, therefore, do not contribute to 
the first-order corrections. We will return to this point shortly. 



When particularized to FLRW cosmologies for matter with constant equation of state P = wp, the time-time 
component of the field equations (6) leads to 
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where H = d/a, and we have set K = for simplicity. From this equation, one readily sees that H 2 vanishes when 
fg — » 0, which corresponds to ^ = —Rp/2 (see [16] for a thorough discussion of the conditions at which can 
vanish). If > 0, this condition is satisfied by sources with w > 1/3. If Rp < 0, then the condition is satisfied by 
w < 1/3. The vanishing of H 2 in this model occurs when a > and implies that the expansion factor has reached a 
minimum. This means that the Universe started in a contracting phase, reached a minimum size, and then bounced to 
a new expanding phase, as shown in Fig. 3. From this plot one can verify that at low curvatures the expansion factor 
follows very accurately the GR solution. It is just very near the characteristic curvature/density scale of the theory 
where the solution departs from GR and the singularity is avoided. The reason for this dynamical change can be found 
in the non-perturbative effects introduced by the terms V v fp and d^/ndy/p on the right hand side of (6), which are 
multiplied by fg l and fj[ 2 , respectively. Though, as we advanced above, the latter contributions do not play any role 
at the perturbative level, they dominate the right hand side of the equations and lead to a strong backreaction that cures 
the singularity when the characteristic curvature/density scale of the theory is reached. We thus see that in ultraviolet 
modified models there exist non-perturbative mechanisms able to protect the theory from certain singularities. Whether 
or not these mechanisms appear in Universes with less symmetry or in other critical situations such as in black hole 
scenarios is a matter which deserves further study. 



In this work we have studied two generic Palatini f(R) models: one modified in the infrared, f(R) = R — fJ. 4 /R, 
and another modified in the ultraviolet, f(R) = R+R 2 /Rp. The infrared corrected model, initially considered as a 
potential candidate to explain the cosmic speedup, has been shown to be ruled out by observations. Models like this 
one are characterized by infrared curvature corrections which cause strong backreaction effects in regions of very low 
matter-energy density such as the zeros of atomic wavefunctions. We have shown that the ground state of Hydrogen 
can be studied using perturbative methods, which guarantees the validity of our results, and turns out to be unstable. 
Other excited states exhibit non-perturbative behavior and confirm that strong modifications of the physics at such 
scales should occur, which supports our conclusions. 

We would like to remark that the 1 /R model in Palatini formalism had not been convincingly ruled out until the 
publication of [12]. This model, first proposed by Vollick [15] in Palatini formalism, was immediately criticized 
by Flanagan [17] arguing that it could be ruled out according to our knowledge on electron-electron scattering 
experiments. Despite Flanagan's efforts to carry out explicit calculations to verify that claim and respond to Vollick's 
comment [18], he concluded [19] that strong coupling effects could make impossible the calculation of observable 
physical effects at experimentally accessible energies. The theory could thus not be ruled out on grounds of electron- 
electron scattering. With the analysis of the Hydrogen atom presented here [12], whose ground state can be treated 
fully perturbatively, we believe that this discussion has been finally settled. 

On the other hand, ultraviolet corrected models manifest interesting features in the context of the very early 
Universe. Non-perturbative mechanisms, analogous to those that spoil infrared corrected models, act in this case but 
with a successful result: the Big Bang singularity can be removed. In the cuadratic model considered here, the solution 
of the modified equations represents a deformation of the general relativistic solution which depends on the parameter 
Rp. (This is so because in Palatini f(R) there are no new degrees of freedom and the equations are of second-order, 
like in GR). If this parameter is set to infinity, then the singular solution of GR is recovered. For any finite value of 
Rp the cosmic singularity is avoided 1 (for some equations of state). Furthermore, the avoidance of the singularity is a 
non-perturbative effect which could not be guessed from a perturbative analysis at low energies. This result is rooted 



1 The quadratic Lagrangian f(R) = R + R 2 /Rp has been criticized in [20] on the basis of the existence of surface singularities in spherically 
symmetric stellar solutions. However, for Planck scale corrected models such singularities could never exist in nature (see [21] for details). 



SUMMARY AND CONCLUSIONS 



in the new role played by the matter sources on the geometry due to its coupling to the metric via the connection. We 
thus see that Palatini f(R) theories provide new gravitational mechanisms that increase the role of the matter sources 
in the spacetime dynamics. 

Though Palatini f(R) theories were originally motivated by phenomenological considerations (cosmic speedup), 
today we know that they may be suitable for describing certain properties of quantum spacetimes. In fact, it was 
recently discovered [22] that the effective equations of loop quantum cosmology [23], an approach to quantum 
cosmology based on non-perturbative quantization techniques, can be derived from a Palatini f(R) theory with high 
curvature corrections. This theory leads to a non-singular bouncing cosmology, and the mechanism responsible for the 
bounce is analogous to that acting in the R + R 2 /Rp model considered here. Our results support the need for further 
analyses of the dynamics of Palatini theories. 
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